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04 . Abstract 

The potential for studying CP violation in neutrino oscillations 
! using conventional and beams is examined. For Am^ << 

CN \ ^ m 3i an d fixed neutrino energy, E u , the CP violating asymmetry 

A = Pfa -> u e ) - -> z/ e )/P(z/ M -> i/ e ) + P(z/ M -> v e ) in vac- 
uum is shown to be measurable with roughly equal maximal sta- 
tistical precision at distances L n ~ (2n+l) ( 2 . 7r£ j / ) , n = 0, 1, 2 . . . 

OO ■ (up to some n in the leading Am^ approximation). For extra 

long baselines, n > 1, the falloff in detected oscillation events, 
O N = N Ve + N Pe , by ~ l/(2n + l) 2 is compensated by a factor 

~ (2n+ 1) increase in the asymmetry such that the statistical fig- 
ure of merit F.O.M. = (5A/A)~ 2 = A 2 N/1 — A 2 is approximately 



(N 



o independent of n. However, for the larger n > 1 asymmetries, 
P^. some backgrounds as well as systematic uncertainties from flux 

normalization, detector acceptance etc. which scale with distance 
as l/(2n + l) 2 are shown to be relatively suppressed in 5 A/ A 
by ~ l/(2n + 1). Also, low energy E u ~ O(lGeV) extra long 
^ ■ baseline experiments with L n ~ 1200-2900 km (corresponding to 

Amjj ~ 3 x 10~ 3 eV 2 and n = 1-3) are better matched to the 
remote locations of some proposed very large proton decay detec- 
tors which would be essential for neutrino CP violation studies. 
Effects of matter and realistic neutrino beam energy spread on 
extra long baseline CP violation experiments are briefly discussed. 

Experimental evidence for neutrino masses, mixing and oscillations has 
become compelling. Atmospheric neutrino flux measurements by the Super- 
Kamiokande (SK) collaboration Jl) indicate large mixing — > v T oscillations 
governed by a neutrino mass squared difference 

Am 2 2 = m\ -m\ ~ (1.6 -4) x 10~ 3 eV 2 (1) 
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(The sign of Am| is not yet determined, but here we assume m 3 > 777.2.) 
Also, recent results from SNO on the solar v e flux when combined with 
earlier SK and other || solar neutrino data strongly point to the so-called 
Large Mixing Angle (LMA) scenario with Q 

Am 2 21 = m 2 - m 2 ~ (2 - 10) x l(T 5 eV 2 (2) 
In terms of the 3x3 mixing matrix 
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Ci = cos6i , Sj = sin^j , 2 = 1,2,3 



those combined results imply relatively large mixing 



sin^ ~ 0.4 (4) 
sin 6 2 ~ 0.707 

On the other hand, failure to observe oscillations in u e — > v x reactor searches 
gives the bound [|] (for Am^ ~ 3 x 10" 3 eV 2 ) 

sin 2 2^ 3 <0.12 or s 3 < 0.18 (90% CL) (5) 

suggesting small first-third generation mixing. In the case of the phase, 5, 
nothing is currently known. Determination of 5 and 83 will tell us the amount 
of CP violation inherent to lepton charged current mixing via the Jarlskog 
CP violation invariant |7J 

J cp = siS2S3Cic 2 c 2 sin5 (6) 

Future experimental efforts will aim to measure the above parameters 
with high precision and search for potential deviations from expectations due 
to "new physics" such as sterile neutrino mixing. Towards those ends, the 
KamLAND reactor experiment || should pin down the values of both Am^ 
and sin#i, if the LMA solar neutrino solution is indeed correct. Similarly, 
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K2K and MINOS should better determine Am| 2 and sin6» 2 via 
disappearance. Further down the road, a recent proposal for a JHF- 
Kamioka — > i/ e oscillation experiment using a conventional horn generated 
pion decay neutrino beam of ~ 0.7 GeV from their future 0.77 megawatt 
proton synchrotron and the large (22kton) SK detector located 295 km away 
would search for sin 2 2# 3 down to 0.006 (a factor of 20 improvement over the 
current bound) and measure Am| 2 to ±10~ 4 eV 2 and sin 2 26> 2 to about ±1%. 

The reach of the JHF-Kamioka proposal || illustrates the significant 
potential of intense low energy conventional neutrino beams for studying 
oscillations when a very large detector such as SK is available. The next 
generation of proton decay detectors will have fiducial volumes 20-40 times 
larger than SK and will probably provide our first opportunity to probe CP 
violation in neutrino oscillations and thereby determine the phase 5. Indeed, 
the JHF-Kamioka neutrino project 0] envisions a phase II upgrade of its 
accelerator complex to 4 Megawatts, which when combined with a possible 
1000 kton Hyper-Kamiokande detector, could increase their yearly statistics 
by a factor of 200. Studies H indicate that with such a facility, in 2 years 
of and 6 years of running, they would statistically probe CP violation 
well below 5 ~ 7r/9, assuming systematic uncertainties in the comparison of 
v n — > v e and P,x —> v e and backgrounds can be controlled. 

In this paper, we examine potential advantages of studying CP violation 
with conventional and beams at extra long baseline distances 

L n ~(2n + 1)L , n = 1,2,3... (7a) 

which are larger by odd multiples 3,5,7 etc. than the generally considered 
first oscillation maximum 

As we will show, to leading order in Am^/Am^ << 1, the statistical figure 
of merit (F.O.M.) 

F.O.M. = {5A/A)- 2 (8) 

which describes ones ability to measure A, Jqp and sin 8 is roughly inde- 
pendent of n (up to some n where the leading Am^ approximation starts 
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to fail). That finding implies that CP violation experiments at the first 
few L n ,n = 0,1,2..., are (roughly) statistically equivalent to first approx- 
imation (until one runs out of events). However, the extra long baselines 
with n > 1 have a potentially significant advantage when it comes to back- 
grounds and systematic uncertainties from neutrino flux normalization, de- 
tector acceptance etc. which scale down by l/(2n + l) 2 with the larger 
distances. Such effects are relatively suppressed in 5 A/ A by ~ l/2n + 1. 
Furthermore, the longer distances may be better matched to the remote 



locations of very large proton decay detectors JLT] which are unlikely to 
be located 413(E l/ /GeV) ( 3X A„J eV2 ) km from existing intense sources 
of E v ~ (9(1 GeV) neutrino beams such as the AGS at Brookhaven or the 
Fermilab Booster. For example Brookhaven- Homestake = 2540km Fermilab- 
Homestake = 1290km while Brookhaven-Carlsbad = 2920km and Fermilab- 



Carlsbad = 1770km ||12| . (An exceptional case is the JHF-Hyper-Kamioka 
proposal || which is fixed at ~ 295 km and thus optimized for and con- 
strained to E u ~ 0.7 ( 3x i°J eV2 ) GeV CP violation studies.) 

We begin by giving the exact — > v e oscillation probability in vacuum 
at a distance L for neutrino energy E v (matter effects are subsequently con- 
sidered). In terms of the mixing angles in ([3D, one finds (for — > v e and 

V e — >■ Vp, JcP — > —Jcp) 

P(vtt v e) = ^(4 s l c l + j cp sin A 2 i) sin 2 4p 

+2(siS2S3CiC2C 2 cos S — s\s\s\c^) sin A 3 i sin A 2 i 

+4(s 1 c 1 c 2 c 3 + s x s 2 SzC z — 2s 1 S2S3CiC2C 3 cos o — Jcp sin A 31 ) sin 

, o/ 2 r 2 2 2 2\ ■ 2 ^31 . 2 ^21 

+o( s i s 2 s 3Ci c 2 c 3 cos o — S1S2S3C3J sin ——sin — ^— 
where Jcp is given in @ and 

A 31 = Am 2 31 L/2E u 

A 21 = Am 2 21 L/2E u (10) 



For studies of CP violation, one considers the asymmetry 



A ~ P(yU>1 — ^ - P(yUfl — ^ fill 



To leading order in A 2 i (assumed to be small), one finds 



P(^ -> v e ) ~ 4^ S 2 C 2 sin 2 + 0(A2i) (12fl) 
A c Jcp ™f 21 * 2SlClC2Si ^ f + 0«) (126) 



The asymmetry grows linearly with distance, but for fixed detector size and 
neutrino energy, the flux of neutrinos decreases as ~ 1/L 2 . The experimental 
determination of A will entail measuring N Ve events from — * v e and N Pe 
from — > z/ e oscillations such that 

AT - r N- 

A "° = t^t, (13) 

where r is the ratio of [y^ flux) x (z/ e cross-section) x (u e N — > e _ iV' acceptance) x 
(running time) divided by the same quantity for antineutrinos. The running 
times would normally (for small asymmetries) be arranged to make r = 1 
up to systematic uncertainties. (For large asymmetries, a different running 
strategy might be employed in which run times are further relatively weighted 

by H-) 

For our statistical discussion, we set r = 1 and find the statistical figure 
of merit (F.O.M.) regarding how well the asymmetry and therefore sin S are 
measured is given by 



F.O.M. = (5A/A)- 2 = A 2 N/1 - A 2 

N = N Ue + Np e (14) 

Note that it is linear in N and quadratic in A. From (12), one finds to leading 
order in A 2 i that the statistical F.O.M. is optimized at the oscillation maxima 

Aqi AmlL . . . . 

~f = ""dlr = (2n + 1)71/2 (15) 
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At distances L n , the growth in A by a factor of (2n + 1) compensates for 
the reduction in N by l/(2n+ l) 2 due to flux reduction, leaving the F.O.M. 
roughly independent of n (at least until the small A 2 i approximation breaks 
down). Where our approximation starts to break down depends on S3 and 
Am^/Arngj as our examples will subsequently illustrate. Note also, that in 
our approximation the F.O.M. for 5A/A is insensitive to the value of S3, at 
least until S3 becomes so small that the approximation in (12a) breaks down. 

To illustrate the above finding and its domain of validity, we consider 
measurements carried out at different L n for E v = 1 GeV beams employing 
the following realistic parameters 

si = 0.4 ci = 0.916 

s 2 = 0.707 c 2 = 0.707 

s 3 = 0.10 c 3 = 0.995 

sin 5 = 0.707 cos S = 0.707 

Am 2 31 = 3 x 10~ 3 eV 2 Am 2 21 = 5 x 10~ 5 eV 2 

in the full expression of eq. @. For definiteness, we assume running times (2 
years of and 6 years of v^) , fluxes and detector assumptions given in phase 
II of the JHF-Hyper-Kamiokande proposal || which envisions a 1000 kton 
water cerenkov detector and 4 megawatt proton synchrotron. We assume 
cuts on the events that reduce the oscillation signals by about a factor of 
2. In that way, our statistics roughly represent simple scaling of the JHF- 
Hyper-Kamiokande proposal. In an actual experiment unconstrained by the 
295 km requirement, one would probably employ somewhat higher (average) 
energy neutrino beams, looser cuts and perhaps a wide band beam at larger 
distances. However, those potential gains in statistics might be somewhat 
offset by oscillation losses at larger distances (see subsequent discussion); 
so, our table numbers should be fairly realistic but are not to be taken too 
literally. 

Table [l] nicely illustrates the main point of this paper. The statistical 
F.O.M. is rather constant for n = 0-3. Beyond those distances, the Am 21 
governed oscillation starts to become significant and eventually begins to re- 
duce the asymmetry (although it remains quite appreciable). Of course, at 
the larger distances, matter effects and beam characteristics become more 
important and should be carefully folded into any distance-beam energy op- 
timization. 
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n 


L n (km) 


K e + N 9e 


A 


F.O.M. 





413 


16000 


0.134 


293 


1 


1239 


1952 


0.366 


302 


2 


2065 


824 


0.516 


299 


3 


2891 


512 


0.586 


268 


4 


3717 


384 


0.601 


217 


5 


4543 


315 


0.585 


164 



Table 1: Number of total oscillation events (after cuts) N Ue + N Pe , the CP 
violating asymmetry A, and F.O.M. at different distances for the mixing 
parameters given above. E v ~ 1 GeV, r = 1 and the experimental conditions 
given in the text are assumed. A F.O.M = 293 corresponds to about a 17 
sigma determination of the asymmetry A. 



To demonstrate the robustness of our result, we give in table || changes 
brought about by reducing S3 from 0.1 to 0.05. Note, that the oscillation 
statistics are reduced, but the F.O.M. is rather unaffected (at least for n = 0- 
2). Here the dilution due to A21 oscillations starts to set in earlier. 



n 


L n (km) 


K e + N Pe 


A 


F.O.M. 





413 


4168 


0.257 


295 


1 


1239 


616 


0.580 


312 


2 


2065 


328 


0.650 


240 


3 


2891 


248 


0.600 


140 


4 


3717 


216 


0.536 


87 



Table 2: Parameters the same as table [l] except S3 = 0.05 instead of 0.10. 

The measurement of CP violation as illustrated in tables [1] and ^| is sta- 
tistically significant, ~ 17 sigma for the first few n. One could half the 
running time, descope the large water cerenkov detector to a more man- 
ageable 500kton, go from a 4 to a 1 Megawatt proton synchrotron and still 
have a robust 4 sigma signal. However, a reduced asymmetry due to a much 
smaller sin^ ^l 1 than the 0.012 value assumed here would require the higher 
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statistics. One should have a better idea regarding experimental demands as 
Am| 1; Arrijj and the 9i are more precisely determined. 

Where the extra long baselines have a potential significant advantage is 
in the area of systematic errors from flux normalization, detector acceptance 
and those backgrounds that scale with distance as l/(2n+ l) 2 . To illustrate 
that point, we consider the experimental asymmetry in ([13]) where r ~ 1 
contains uncertainties from flux, cross-section, acceptance and other nor- 
malization effects. The fractional systematic error 5r/r translates into the 
following asymmetry error 

5A/A ~ - — — - (16) 

For larger asymmetries, the effect of Sr/r is significantly reduced. For asym- 
metries exhibiting the (2n + 1) growth, such systematics are reduced by 
~ l/2n+ 1. This advantage can be significant, particularly if the CP violat- 
ing asymmetry turns out to be much smaller than in our previous illustrative 
examples. 

Other backgrounds, such as v e beam contamination and misidentified 
neutral current residual events remaining after cuts, that scale down by 1/L 2 
will similarly be reduced in relative importance by ~ 1/A for larger asym- 
metries. Non-distance dependent effects such as cosmic ray backgrounds 
would of course be worse for the reduced statistics of extra long baseline ex- 
periments. Fortunately, they are generally very small in deep underground 
experiments because of directionality and the precise timing constraints im- 
posed by beam pulsing. Nevertheless, such effects are potentially more dan- 
gerous and must be carefully studied. 

For realistic experimental considerations, the full oscillation formula in 



(H) supplemented with matter effects |L3| must be employed. In addition, 
neutrino beam energy spread, detector acceptance, backgrounds etc. must 
be considered. Here, we make some brief qualitative observations concerning 
such effects. They will be further examined in a subsequent more detailed 
study. 

At extra long distances, L n , n = 1-3, the smaller frequency A 2 i os- 
cillations and matter effects will effectively shift somewhat the oscillation 
peaks. More important, matter effects will induce a fairly significant CP 
non-violating (positive) asymmetry which must be well controlled and un- 
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derstood in any serious study. For the distances considered here ~ 1200-2900 
km and neutrino energies E v ~ 0(1 GeV), such calculations can be carried 
out in a straightforward way that exactly includes Am^, Am^ and matter 



density effects [[L4| . Such a study has been nicely illustrated in ref. 15. There, 
one finds for a distance of 2900 km that matter effects are most important 
for the n = peak at E v ~ 6 GeV. At lower energies, v ^ oscillation peaks 
are also somewhat shifted but less amplified due to Am^ and matter effects. 
Those shifts may influence the experimental beam energy setting for a given 
baseline distance, but they should not significantly change our discussion 
about the statistical significance of extra long oscillation studies near the 
first few L n with n = 1, 2, 3 as long as E v is not too large. The CP violating 
asymmetry increase at the longer distances should continue to compensate 
for the flux falloff in the F.O.M. at least for the first few n > 1. If the matter 
asymmetry and CP violating asymmetry have the same sign (both positive), 
we can have potentially enormous combined asymmetries at long distances 
with the I— A 2 in the F.O.M. denominator playing a significant role. Alterna- 
tively, for opposite signs, matter and CP violation asymmetries will partially 
cancel. Of course, such a cancellation is also possible and potentially more 
problematic for n = 0. (In the JHF-Kamioka discussion of ref. 5, a partial 
cancellation of matter and CP violating asymmetries which dilutes the sig- 
nal is apparent for the phase chosen as an illustration 5 = —45°.) Generally, 
one expects better control of matter induced asymmetries for the larger CP 
violating asymmetries of extra long, n > 1, baselines. 

Also of importance for extra long beamline experiments is the neutrino 
energy spectrum [|l(|. For conventional pion decay generated neutrino beams 
from 16 ~ 28 GeV proton drivers, one can either employ a narrow band 
(energy spread) beam centered at some energy E v ~ 1 GeV by placing the 
detector off-beam axis by some small angle or opt for a more intense wide 
band beam with considerable support in the E v ~ 0.5-2 GeV range |5]. 16 



Because the oscillation peaks at n > 1 have a narrower energy width, it will 
be more challenging to tailor a narrow band beam to the detector distance 
and oscillation parameters as n grows. Instead, the wide band neutrino 
beam may have advantages. In addition to the higher flux of neutrinos, the 
energy spread would simultaneously cover several oscillation maxima. For 
example, if the E v ~ 2 GeV component corresponds to the n = 1 peak 
(for Li ~ 2500km), there will also be potentially observable peaks at E v 
approximately 1.2 GeV (n = 2), 0.86 GeV (n = 3) and 0.67 GeV (n = 4). Of 



9 



course, there will also be oscillation minima between the peaks that would 
reduce statistics (roughly wiping out some flux gains of a wide band beam). 
Nevertheless, an experiment that covers and can partially unravel the n = 1- 
4 oscillation peaks has the potential (if statistics suffice) to not only probe CP 
violation, but to determine the sign of Am^ and study the changing matter 
effects in detail. Such potential capabilities provide strong motivation for 
extra long baseline experiments. Those features will be explored in detail in 
a subsequent study. 



References 



[i] 

[2] 
[3 



[4] 

[5] 

[6] 

[7] 
[8 

[9 
[10 



Super-Kamiokande collaboration, Phys. Rev. Lett. 81, 1562 (1998); 

Q.R. Ahmad et al, (SNO Collaboration), Phys. Rev. Lett. 

B.T. Cleveland et al, Astro. J. 496, 505 (1998); K. Hirata et al, Phys. 
Rev. D44, 2241 (1991); W. Hampel et al, Phys. Lett. B447, 127 (1999); 
J. Abdurashitov et al, Phys. Rev. C60, 055801 (1999). 

J. Bahcall, M. Gonzalez- Garcia and C. Pena-Garay, hep-ph/106258 and 
references therein. 



Y. Itow et al, "The JHF-Kamioka Neutrino Project" , |K"ep-ex/010601~9 



The sign of the Jqp terms in this paper seems to be the opposite of ours. 

CHOOZ: M. Apollonio et al, Phys. Lett. B466, 415 (1999); Palo Verde: 
F. Boehm et al, Nucl. Phys. Proc. Suppl. 91, 91 (2001). 

C. Jarlskog, Z. Phys. C29, 491 (1985); Phys. Rev. D35, 1685, (1987). 

A. Piepke (KamLAND collaboration), Nucl. Phys. B91 (Proc. Suppl.) 
(2000). 



S. Ahn et al, frep-ex/010300T 



MINOS home page |http:// www . hep . anl . gov /ndk / hyp er t ext / minos| . 
htm/ 



[11] C.K. Jung, |hep-ex/0005046 



10 



[12] V. Barger, D. Marfatia and K. Whisnant, |hep-ph/0108090| . 

[13] L. Wolfenstein, Phys. Rev. D17, 2369 (1978); S.P. Mikheyev and A. 
Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 913 (1986)]. 

[14] J. Arafune, M. Koike and J. Sato, Phys. Rev. D56, 3093 (1997); H. 
Zaglauer and K. Schwarzer, Z. Phys. C40, 273 (1988); T. Ohlsson and 
H. Snellman, J. Math. Phys. 41, 2768 (2000). 

[15] I. Mocioiu and R. Shrock, [hep-ph/0106139| . This paper contains many 
references to studies of matter effects. 

[16] D. Beavis et al, BNL Report 52459. 



11 



